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1.  INTRODUCTION 

Let  Xj,  Xg,  ...  be  Independent  identically  distributed  (l.l.d.)  random 
variables  (r.v.'s)  denoting  llfelengths  with  a  common  continuous  distribu¬ 
tion  function  (d.f. )  F,  and  let  Yj,  Yg,  ...  be  l.l.d.  right-censoring  r.v.'s 
with  a  cooinon  continuous  d.f.  H.  We  assume  that  (Xj,  Yj),  (Xg,  Yg),  ... 
is  an  l.l.d.  sequence  of  random  pairs  with  Independent  nonnegative  components 
defined  on  a  common  probability  space. 

Let  1(A)  denote  the  indicator  function  of  the  set  A,  let  *  mlnlX^,  Y^> 
and  »  I(X^  <  Y^),  1  »  1,  ....  n.  Ue  consider  the  problem  of  estimating 
Fi  1  -  F,  the  underlying  survival  function,  from  the  sample  (Z^,  4j),  ....  (Zn,  4fl). 

Kltchln,  Langberg,  and  Proschan  (KLP)( 1980a)  introduce  a  new  estimator 
of  F|  the  Piecewise  Exponential  Estimator  (PEXE),  a  competitor  to  the  tradi¬ 
tional  Product  Limit  Estimator  (PLE),  Introduced  by  Kaplan  and  Meier  (1958). 

KLP  show  that  the  PEXE  Is  a  strongly  consistent  estimator  of  F  and  that  the 
standarized  PEXE  process  converges  weakly  to  a  Gaussian  process.  This  Gaussian 
process  has  the  same  covariance  structure  as  that  of  the  standarized  PLE 
process.  [For  details  see  KLP  (1980a)  and  (1980b).]  Thus  the  PEXE  has 
the  same  asymptotic  behavior  as  the  PLE.  [For  the  strong  consistency  of 
the  PLE  see  Peterson  (1977),  Langberg,  Proschan,  and  Qulnzl  (1980)  or  KLP 
(1980a).  For  the  weak  convergence  of  the  standarized  PLE  see  Breslow  and 
Crowley  (1974)  or  KLP  (1980b).]  In  contrast  to  the  PLE,  which  Is  a  step 
function  with  jumps  at  the  observed  failures,  the  PEXE  Is  a  continuous  and 
strictly  decreasing  function  up  to  the  last  failure.  In  many  life  testing 
situations  the  survival  function  Is  anticipated  to  be  strictly  decreasing 
smoothly  over  time.  Thus,  In  many  situations  the  PEXE  will  be  the  more 
appropriate  than  the  PLE.  Further,  the  value  of  the  PEXE  In  an  Interval 


between  any  two  successive  failures  depends  on  the  actual  withdrawal 
times  up  to  the  failure  time  that  determines  the  right-hand  side  of  the 
Interval,  and  on  that  failure  time.  In  contrast  the  value  of  the  PLE 
between  two  successive  failures  depends  on  the  number  of  withdrawals  up 
to  the  failure  that  determines  the  left-hand  side  of  the  Interval,  and 
does  not  depend  on  the  right-hand  side  failure.  [For  further  details 
see  KLP  (1980a).] 

In  this  paper  we  consider  the  small  sample  properties  of  the  two  esti¬ 
mators.  Our  calculations  are  performed  under  the  assumption  that  the  life- 
length  X  and  tlie  time  to  censorship  Y  have  proportional  hazards.  (See  Deflnl 
tlon  2.3.) 

In  Section  2  we  define  the  PEXE,  the  PLE  and  the  proportional  hazards 
models.  Then  we  present  a  useful  characterization  of  these  models  that  Is 
used  In  the  last  two  sections. 

In  Section  3  we  present  formulas  for  the  mean  and  variance  of  the 
PLE,  and  upper  and  lower  bounds  for  the  mean  and  variance  of  the  PEXE. 

All  computations  In  Section  3  are  conducted  under  the  proportional 
hazards  assumption.  In  Section  4  we  compute  the  mean  and  variance  of 
the  PLE,  and  corresponding  upper  and  lower  bounds  of  the  PEXE.  We 
compare  our  exact  result  for  the  bias  of  the  PLE  with  Efron's  (1967) 
upper  bound.  It  Is  seen  that  Efron's  bound.  In  the  proportional  hazards 
model.  Is  quite  loose.  We  also  compare  tlie  exact  variance  of  the  PLE 
with  an  approximate  variance  given  by  Kaplan  and  Meier  (1958).  The 
latter  approximation  Is  found  to  be  quite  good. 


It  Is  also  seen  In  Section  4  that  our  upper  bound  for  the  mean  of  the 
PEXE  Is  reasonably  close  to  the  true  value  F(t)  but  the  lower  bound  Is  not 
good.  Furthermore,  the  upper  bound  for  the  variance  of  the  PEXE  Is  compared 
to  the  approximation  given  by  Kaplan  and  Meier  (1958)  for  the  variance  of 
the  PLE.  (Recall  that  the  PEXE  and  the  PLE  have  the  same  asymptotic  vari¬ 
ance,  and  thus  the  Kaplan-Meler  approximation  can  also  be  viewed  as  an 
approximation  to  the  (small -sample)  variance  of  the  PEXE.)  The  upper 
bound  Is  seen  to  be  too  conservative,  and  tighter  bounds  are  needed. 


2.  PRELIMINARIES 


In  this  section  we  define  the  PEXE,  the  PLE,  and  the  proportional  hazards 
models.  Then  we  present  a  characterization  of  these  models. 

First  we  Introduce  some  notation.  Let  ^(t)  ■  >  t),  t€[0,  •), 

be  the  nunber  of  observations  In  the  sample  at  time  t,  and  let  T(n)  *  E?«l*1 
be  the  total  number  of  failures.  Further,  let  ...  <  w(T(n) )  denote 

the  consecutive  observed  failures  with  =  0,  and  for  t(n)  >,  1,  let 

i »  _ 

r4  *  [/  '  '  nK-fuJdu]"1  be  the  number  of  failures  per  unit  time  in  the  interval 

i  ii  n 

W(1-l) 

j]»  i  *  1»  »•»»  T(n). 

He  are  ready  to  present  the  definitions  of  the  PEXE  and  the  PLE. 

Definition  2.1.  For  T(n)  >  1  let  a{  *  (ll^j  -  ))r^ ,  1  ■  1 . T(n). 

Then  the  Piecewise  Exponential  Estimator.  (PEXE),  denoted  by  5n(t),  Is  equal 
to  1  on  the  set  (T(n)  *  0  or  te(-«,  0]},  Is  equal  to  expl-rj”^  -  (t  -  w(i_i))r(i) 
on  the  set  (x(n)  £  1,  t«(U^_j  j,  W^j]},  1*1,  ...,  x(n),  and  Is  equal 
to  exp{-£jij^>,  on  the  set  (x(n)  >  1,  t€(M(T(n))»  •))• 

Definition  2.2.  Let  ?n(t»)  *  n-1£?sjI(Z^  >.  t),  t<[0,»).  Then  the 

Product  Limit  Estimator.  (PLE),  denoted  by  Ffl(t),  is  equal  to  1  on  the  set 

M  . 

(x(n)  *  0,  or  tc(-«,0)}.  Is  equal  to  n  [Kn(lJ^j  j)][Kn(H^  j-)]”1,  on  the 

set  (x(n)  >_  1,  t«[U/|_jj,  W^j)},  1*1,  ...,  x(n),  and  Is  equal  to 
»(n)_  -  "  -1 

»  Cl^>{M(J))3CKn(W{J)-)]  \  on  the  set  <x(n)  >  1,  tc[W(T(n)),  -)}. 

He  note  that  originally  Kaplan  and  Meier  (1958)  left  the  PLE  undetermined 

on  the  set  (x(n)  >  1,  t«[  max  Z.,«)).  (See  also  Peterson  (1977).) 

l<1<n  1 

Next  we  define  a  proportional  hazards  model.  Let  (Uj,  Ug)  be  a  pair  of 
Independent  nonnegative  r.v.'s  with  continuous  d.f.'s  Gj,  Gg  respectively. 

Further,  let  ^  *  1  -  G^,  •  sup(t:  Gq(t)  <  1),  and  let  R^(t)  -  -tri^(t), 

t«(-»,  a(G  )),  be  the  hazard  function  corresponding  to  the  d.f.  G ,  q  •  1,  2. 
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Definition  2.3.  We  say  that  the  pair  (Uj,  U^)  follows  a  proportional 
hazards  model  If  there  is  a  positive  real  number  0  such  that: 

£>(t)  *  tGjtt )3B  for  t«C0,-).  (2.1) 

Note  that  from  (2.1)  a(G^)  *  a(Gg),  and  that  the  hazard  functions  corresponding 
to  Gj,  Gg  are  proportional : 

Rx(t)  -  8R2(t),  t«[0,  a(Gj )). 


Next  we  characterize  the  proportional  hazards  models. 

Theorem  2.4.  The  pair  (Uj,  U2),  0  <  PUJj  <.  U2>  <  1,  follows  a  propor¬ 
tional  hazards  model  If  and  only  If 

The  r.v.'s  U  =  m1r.{Uj,  Ug),  and  €  *  I(Uj  ,<  U2)  are  independent.  (2.2) 

Proof.  Let  G  be  a  continuous  d.f.,  let  E  =  1  -  G,  and  assune  that 
G(0)  *  0.  Then  for  e«[>l,«),  and  ze[0,a(G)): 


J[S(u)]°dG(u) 

2 


kl  +  oylCG(z)l*+lt  e  >  -l 


£ 


(2.3) 


-ih£(z)  ,  e  *  -l. 

First  we  prove  that  (2.1)  Implies  (2.2).  To  verify  (2.2)  It  suffices 
to  show  that  for  z«C0,») 


P(U  >  z,  C  ■  1)  «  P(U  >  z)PU  -  1).  (2.4) 

By  (2.1)  and  (2.3) 

P(U  >  z,  C  ■  1)  *  /  £g(u)dGj(u)  «/  CGj(u)]*dGj(u) 

•  (1  ♦  sr^fl^U)]0*1  ■  (1  ♦  B^PW  >  z). 
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In  particular  P{£  «  1}  *  P{U  >  0,  e  ■  1}  »  (1  ♦  p)”*.  Thus,  (2.4)  holds. 
Consequently  (2.1)  Implies  (2.2). 

Now  we  show  that  (2.2)  implies  (2.1).  Let  y  *  P{Uj  <  Ug).  By  (2.2): 

Y^fz^z)  •  /^(uJdGjfu),  z«[0,«).  (2.6) 

Upon  integration  by  parts 

Y/<a^(u)dG2(u)  *  (1  -  Y  If  (^(uJdGj^u),  z«[0,«).  (2.7) 

How  (2.7)  implies  for  z«[0,  m1n{o(Gj),  a(G2))),  that: 

y/C^uH'Wu)  =  (1  -Y  jf[G|(u)]”1dG.(u).  (2.8) 

0  *  c  0  1 

By  (2.3)  and  (2.8),  for  tcCO^InMG^.afGg))), 

-  CE1(t)f'1(1'r).  (2.9) 

Consequently  (2.2)  implies  (2.1).  || 

Since  the  random  pairs  (X^,  )»  1  *  1,  2,  ...  are  l.l.d.  with  Inde¬ 

pendent  components  it  follows  that: 

The  r.v.’s  1  ■  1,  2,  ...  are  l.l.d.  Bernoulli  with  the 

1  (2.10) 

parameter  y  ■  P(Xj  <_  Y^). 

Further,  under  the  assumption  that  (Xj,  Yj)  follows  a  proportional  hazards 
model  vie  conclude  from  Theorem  2.4  that: 

The  sequences  of  r.v.  s  6. ,  *•*,  and  2.,  Z.,  ... 

14  1  (2.11) 

are  Independent. 

Thus,  the  proportional  hazards  assumption  yields  a  nice  structure  on  the 
random  pairs  (Xj,  Yj),  (X2,  Y2),  ...  .  This  structure  plays  an  Important 
role  In  our  computations,  presented  In  Sections  3  and  4.  Without  this 
assumption  finite  sample  computations  are  extremely  difficult.  (Thus 
our  assumption  Is  also  pragmatic.) 
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3.  MOMENT  CALCULATIONS  AND  BOUNDS  UNDER  PROPORTIONAL  HAZARDS. 

In  this  section  we  present  formulas  for  the  mean  and  variance  of  the 
PLE.  Then  we  give  upper  and  lower  bounds  for  the  mean  of  the  PEXE  and  an 
upper  bound  for  the  variance  of  the  PEXE.  These  are  derived  under  the 
assumption  that  (Xj,  Yj)  follows  a  proportional  hazards  model.  *  Throughout 
we  define  a  sum  and  a  product  over  an  empty  set  of  Indices  as  0  and  1 
respectively. 

First  we  present  formulas  for  the  mean  and  variance  of  the  PLE.  Let 
Kn(t)  -  1  -  K^t),  Tc(t )  -  F(t  )TT(t),  and  K(t)  -  1  -  F(t),  tc[0,«).  Note 
that  by  the  continuity  of  F  the  formula  for  the  PLE,  given  In  Definition 
2.2,  reduces  to: 

V*)  6 

F„(t)  -  n  C(n  -  1  )(n  -  1  +  l)-1]  \  t€(-,.).  (3.1) 

n  1-1 

Let  t«[0,-).  Then  KR(t)  is  a  binomial  r.v.  with  parameters  n  and  K(t), 

61*  ***  5n  are  *•*•<*•  Bernoulli  r.v.'s  with  parameter  y,  and  by  (2.11)  Kn(t) 
Is  Independent  of  fij . <n.  Thus,  by  (3.1),  for  t,  «€((),•): 

Kn(t) 

EfT„(t)>a  «  E  fl  [Y(n  -  1  )a(n  -  1  +  1)"°  +  1  -  y] 
n  1-1 

-  E^(;)CK(t)m(t)3n-qin  (Y(n  -  1)a(n  -  1  +  l)**  +  1  -  y). 
Consequently  for  t«[0,«): 

Efn(t)  -  E5.0(J)[K(t)]qU(t)]n“^Ml  -  Y(n  -  1  ♦  I)-1).  (3.2) 

Var(F  (t)>  -  i"  0(;)[K(t)]q(J(t)3n"q  n  (1  -  Y(2n  -  21  +  l)(n  -  1  ♦  1)“2> 
“  U  “  1*1 

,  (3-3) 

(i  -  t(»  -  i  *  D'1)!2. 


and 
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Equations  (3.2)  and  (3.3)  are  conveniently  used  for  numerical  comparisons 
In  Section  4. 

How  we  obtain  upper  and  lower  bounds  for  the  mean  of  the  PEXE  and  an  upper 
bound  for  the  variance  of  the  PEXE.  note  that 

e"x  1  <  x(l  +  x)"1  <  e"*x+1)  \  X€(0,«),  (3.4) 

and  that  on  the  set  {r(n)  >  1} 

CnKjj(ll(i_i ) )]  £  4^  [nK^(U(i))3  ,  1  =  1.  •••»  x(n).  (3.5) 

For  te(0,»),  let  a(n,  t)  »  max(q:  q  =  0,  ....  x(n),  j  <  t)  be  the  Index 
of  the  largest  preceding  t,  and  on  the  set  {a(n,  t)  <  x(n)>  let 

b(n,  t)  *  m1n(q:  q  *  ?n(t)  ♦  1,  ....  n,  *  1}  be  the  Index  of  the  smallest 
observation  that  failed  following  t.  By  the  formula  for  the  PEXE,  given 
In  Definition  (2.1),  and  by  (3.4),  (3.5): 

3n(t)  <  exp{-E^J»t^i) 

<  exp{-l^J»t)[n?(|(W(i-1))r1> 

ia(AWu(M)W^(M(i-i),  +  n’1^ 


and 


Kn(t)  « 

n  {(n  -  1  ♦  l)(n  -  1  ♦  2)*1)  *  for  t«(0,*), 
1-1 


5„(t)  >  Cl(a(n,  t)  -  v(n))  ♦  I(a(n,  t)  <  T(n))exp{-Ab^n>t j)]exp{-E*ij,t^i) 
>[!(«(",  t)  -  t(n))  ♦  I(a(n,  t)  <  T(n))exp{-ab^n  t jHexpl-E^J^CnK^W^  j)]*1) 
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1  [I(a(n,  t)  *  t(n))  +  |rl(a(n,  t)  <  t(n),  b(n,  t)  <  n  -  1] 

K„(t)  . 

I(Kn(t)  <  n)  nn  {(n  -  i  -  l)(n  -  I)"1)  1 
n  i=l 

i  K  (t )  .  «. 

>  £l(Kn(t)  <  n)  nn^  {(n  -  1  -  l)(n  -  If1}  1  for  t«[0,  -). 

Thus,  for  t,  o«(0,») 

Kn(t) 

E{S„(t)>a  <  E  n  CY(n  -  i  +  1  )°(n  -  1  +  2fa  +  1  -  y3 
n  1*1 


(3.6) 


*  -  1  ♦  2>'“  ♦  »  -  O. 

E{Tn(t)}°  >  r°E  n  [Y(n  -  1  -  1  )a(n  -  if*  +  1  -  y3 


2'aES*J(!)CK(t)3qCK(t)]n'q  !  Cy(h  -  i  -  1  )a(n  -  1)'°  +  1  -  y3. 

9  u  9  -f=i 


Consequently 


#n(t)  <  C nCK(t)]q[K(t)]n-q  I  (1  -  y(«  -  1  ♦  2 r1)  (3.7) 


E\(t)  >  2-1Ej:J[K(t)]<»[K(t)]n-q_Ml  -  y(n  -  i  f1)  (3.8) 


Varff^t)}  -  zqs0(q)CK(t)3qC^(t)]n"^n^(l  -  Y(2n  -21+3) 

(n  - 1  +  2r2)-(2“1i;;J(;)[K(t)3qCK(t)3n_q  2  a  -  Y(n  -  ir1))2. 

9  u  9  1«1 


(3.9) 
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4.  COMPARISON  OF  EXACT  VALUES  WITH  BOUNDS  AND  APPROXIMATE  VALUES 
In  this  section  we  take  the  llfelength  d.f.  F  to  be  exponential  with 
scale  parameter  1  and  the  right-censorship  d.f.  H  to  be  exponential  with 
scale  parameter  8 (p  >  0). 

Table  1  displays  results  obtained  by  using  (3.2)  and  (3.3)  to  obtain 
nunerlcal  values  for  the  mean  and  variance  of  the  PLE.  Then  we  compare 
the  bias  of  the  PLE  to  the  general  bound  for  the  bias  given  by  Efron  (1967), 
namely: 

0  <  E^n(t)  -  F(t)  <  F(t)exp(-nl(t)},  t€[0,»).  (4.1) 

Ue  also  compare  the  exact  variance  of  the  PLE  to  the  approximate  variance 
given  by  Kaplan  and  Meier  (1958): 

Var(F(t)}  -  n-^ttf/CKtu^ujr^Ffu),  t€[0,«).  (4.2) 

0 

We  note  that  [F(t)]^/[K(u)F(u)]"*dF(u)  is  the  variance  of  G(t),  the  Gaussian 

0 

process  obtained  as  the  weak  limit  as  n  ♦  •  of  n1^(Fn(t)  -  F(t)J.  G(t)  is 
also  the  weak  limit  of  n^(5_(t)  -  F(t)l.  [See  KLP  (1980b).]  Thus  the 
right-hand  side  of  (4.2)  can  also  be  considered  as  an  approximation  to 
Var{5n(t)>. 

Ue  also  use  (3.7),  (3.8),  and  (3.9)  to  obtain  lower  and  upper  bounds  fbr 
the  mean  of  the  PEXE,  and  upper  bounds  for  the  variance  of  the  PEXE.  See 
Table  2. 

All  our  calculations  are  performed  under  a  particular  proportional 
hazards  model,  when  X  and  Y  are  exponential  r.v.'s.  However  these  calcula¬ 
tions  apply  to  general  proportional  hazards  models,  provided  the  appropriate 
adjustments  are  made.  More  specifically  assume 
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IT  ■  for  some  0  >  0,  (4.3) 

and  let  R”l(t)  *  1nf{z:  zc[0,a(F)),  -triF(z)  >  t),  t«[0,<»).  Since  F  Is  con¬ 
tinuous  R(R"*(t))  *  t  for  t«[0,»).  Thus: 

The  r.v.'s  R(X)  and  R(Y)  are  exponentially  distributed  with 
scale  parameters  1  and  a  respectively. 

Let  Zi>R  *  mintt^),  R(Y1 ))  and  R  *  I(R(X1 )  <  R(Y- )),  1  »  1 . n.  Then: 

The  random  vectors  <  <Z(1)R,  «(1)>R).  ....  (Z(n)>R,  «(„)>R)  > 

and  <  (R(Z^j),  «^j),  ...»  (R(Z(n))»  >  are  stochastically  equal. 

Wow  let  7n  R  and  R  be  the  PLE  and  PEXE  computed  from  R,  R),  ...» 

^Z(n),R*  6(n),R^»  ■'’••poetlwly.  Then  for  n  «  1,  2,  ...: 

The  two  PLE's  processes  (F_  D(R(t))»  t«[0,»)}  and 

n  |K 

(FR(t),  t<0»«*)>  are  stochastically  equal. 

In  this  section  we  compare  the  exact  bias  and  exact  variance  of  the 
PLE,  for  proportional  hazards  models,  with  Efron's  (1967)  upper  bound 
for  the  bias  and  with  Kaplan  and  Meier's  (1956)  approximate  variance.  Thus, 
we  may  without  loss  of  generality  replace  the  model,  given  by  (4.3),  with 
the  particular  exponential  model  described  In  the  first  paragraph  of  this 
section. 

lie  note  that  the  PEXE  processes  (5„  „(R(t)),  t<[0,«)}  and  (5_(t),  t«[0,-)) 

n*K  n 

are  not  stochastically  equal.  However,  since  the  upper  and  lower  bounds 
of  ?n,  given  In  Section  3,  depend  only  on  Kn(t)  and  T(n),  our  comparisons 
for  the  PEXE  conducted  for  an  exponential  model  apply  to  a  general  proportional 
hazards  model  provided  t  Is  replaced  by  R(t). 
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Next  we  summarize  our  results  for  the  PLE  in  the  following  table. 


1.  Exact  and  Approximate  Values  of  the  Mean  and  Variance  of  the  PLE 


sample  size  n 


mean  of  the  PLE 
bias  of  the  PLE 
Efron's  bound  for  the  bia: 
variance  of  the  PLE 
approximate  variance _ 


=  1.  t 


sample  size  n 


mean  of  the  PLE 
bias  of  the  PLE 
Efron's  bound  for  the  bias 
variance  of  the  PLE 
approximate  variance 


.6065 


.0003 


.6065 


.0107 

0105 


sample  size  n 


mean  of  the  PLE 
bias  of  the  PLE 
Efron's  bound  for  the  bias 
variance  of  the  PLE 
approximate  variance 


.0042 

.0129 

0157 


.0005 

.0107 

0105 


sample  size  n 


variance  of  the  PLE 

approximate  variance 

.0474 

.0423 

,0318 

.0288 

.0236 

.0216 

.0173 

.0186 

.0144 

.0154 

L 

f-2. ±JL  L 

,  ft t)  «  .3679 

sample  size  n 

10 

15 

20 

25 

30 

mean  of  the  PLE 

.4164 

.3946 

.3836 

.3799 

.3744 

bias  of  the  PLE 

.0485 

.0268 

.0160 

.0100 

.0065 

Efron's  bound  for  the  bias 

.3842 

.2995 

.2335 

.1821 

.1419 

variance  of  the  PLE 

.0779 

.0580 

.0462 

.0382 

.0324 

approximate  variance 

.0861 

.0574 

.0344 

.0287 
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We  see  that  Efron's  (1967)  bound  for  the  bias  Is  much  larger  than  the  exact 
bias,  computed  for  proportional  hazards  models.  However  the  approximate 
variance  given  by  Kaplan  and  Meier  (1958)  Is  close  to  the  exact  variance, 
even  for  the  smallest  sample  size  considered  (n  *  10). 

The  bounds  for  the  PEXE  In  Table  2  are  based  on  formulas  (3.7),  (3.8) 
and  (3.9).  The  upper  bound  for  the  mean  of  the  PEXE  Is  reasonably  close 
to  the  true  value  of  F(t)  but  the  lower  bound  for  the  mean  Is  not  close. 

Sharper  bounds  need  to  be  developed. 

For  completeness,  we  Include  the  upper  bound  for  the  variance  of  the 
PEXE  but  until  an  exact  expression  is  obtained  for  the  (small -sample)  vari¬ 
ance  of  the  PEXE,  the  potential  value  and  use  of  this  "pper  bound  Is  diminished. 
However,  vie  can  compare  the  upper  bound  to  the  approximate  variance  given 
by  the  right-hand  side  of  (4.2).  Such  a  comparison  Indicates  that  the  upper 
bound  Is  not  that  good. 

2.  PEXE  Bounds:  Upper  and  Lower  Bounds  for  the 
Mean  and  Upper  Bound  for  the  Variance 


6  -  0.5.  t  «  0.5.  F(t)  «  .6065 


sample  size  n 


10  15 


20 


upper  bound  for  the  mean  of  the  PEXE  .6476  .6348  .6280 
lower  bound  for  the  mean  of  the  PEXE  .2768  .2866  .2911 
upper  bound  for  the  variance  of  the  PEXE  .3649  .3366  .3220 
approximate  variance _  .0274  .0183  .0137 


25  30 

.6239  .6211 
.2937  .2954 
.3131  .3070 
.0110  .0091 


ixlmate  variance 


*  1,  t  «  0.5,  F(t)=  .6065 


sample  size  n 


upper  bound  for  the  mean  of  the  PEXE 
lower  bound  for  the  mean  of  the  PEXE 
upper  bound  for  the  variance  of  the  F 
xlmate  variance  _ 


10 

15 

20 

25 

■a 

6539 

.6391 

.6313 

.6266 

.6233 

2706 

.2831 

.2888 

.2919 

.2940 

3792 

.3463 

.3292 

.3188 

.3119 

0316 

.0211 

.0158 

.0126 

.0105 

sample  size  n 


upper  bound  for  the  mean  of  the  PEXE 
lower  bound  for  the  mean  of  the  PEXE 
upper  bound  for  the  variance  of  the  F 
approximate  variance 


10 

15 

20 

25 

30 

6700 

.6505 

.6400 

.6336 

.6292 

2510 

.2712 

.2810 

.2864 

.2898 

416? 

.3723 

.3488 

.3334 

.3247 

0427 

.0285 

.0213 

.0171 

.0142 

sample  size  n 


=  3679 


10  15  20  25 


upper  bound  for  the  mean  of  the  PEXE  .4453  .4212  .4085  .4007  .3954 
lower  bound  for  the  mean  of  the  PEXE  .1352  .1510  .1597  .1650  .1685 
upper  bound  for  the  variance  of  the  PEXE  .2044  .1723  .1552  .1447  .1376 
approximate  variance _ .0314  .0209  .0157  .0126  .0105 

8  «  1.  t  «  1.  f(t)  «  .3679 


sample  size  n  10  15 


upper  bound  for  the  mean  of  the  PEXE  .4722  .4407  .4237  .4130  .4058 

lower  bound  for  the  mean  of  the  PEXE  .1192  .1365  .1477  .1553  .1606 

upper  bound  for  the  variance  of  the  PEXE  .2386  .1980  .1755  .1613  .1516 

xlmate  variance  .0432  .0288  .0216  .0186 
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8  »  2.  t  ■  1.  F(t)  =  .3679 


sample  size  n 

10 

15  20  25 

30 

upper  bound  for  the  mean  of  the  PEXE 

.5396 

.4966  .4707  .4532 

.4407 

louer  bound  for  the  mean  of  the  PEXE 

.0783 

.0981  .1097  .1205 

.1292 

upper  bound  for  the  variance  of  the  PEXE 

.3247 

.2694  .2364  .2142 

.1282 

approximate  variance 

.0861 

.0574  .0430  .0344 

.0287 

@_=  0,5,  t  =  2,FftJ  «  . 

1353 

sample  size  n 

10 

15  20  25 

30 

upper  bound  for  the  mean  of  the  PEXE 

.2782 

.2375  .2148  .2003 

.1902 

lower  bound  for  the  mean  of  the  PEXE 

.0251 

.0298  .0339  .0374 

.0404 

upper  bound  for  the  variance  of  the  PEXE 

.0938 

.0676  .0544  .0465 

.0412 

approximate  variance 

.0223 

.0155  .0117  .0093 

.0078 

8  - 1. 1  • mi 

«  .1353 

sample  size  n 

10 

15  20  25 

30 

upper  bound  for  the  mean  of  the  PEXE 

.3566 

.3049  .2737  .2525 

.2370 

lower  bound  for  the  mean  of  the  PEXE 

.0170 

.0205  .0234  .0260 

.0283 

upper  bound  for  the  variance  of  the  PEXE 

.1538 

.1123  .0901  .0762 

.0667 

approximate  variance  .0491  .0327  .0245  .0196  .0164 


sample  size  n 

10  15 

20 

25  30 

upper  bound  for  the  mean  of  the  PEXE  .4955  .4398  .4033  .3769  .3565 
lower  bound  for  the  mean  of  the  PEXE  .0043  .0056  .0067  .0078  .0067 
upper  bound  for  the  variance  of  the  PEXE  .2869  .2272  .1915  .1674  .1499 
approximate  variance  _  _ .2457  .1638  .1228  .0983  .0819 
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